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Abstract 

We study the asymptotic behavior of solution of semi-linear PDEs. Neither period- 
icity nor ergodicity will be assumed. In return, we assume that the coefficients admit a 
limit in Cesaro sense. In such a case, the averaged coefficients could be discontinuous. 
We use probabilistic approach based on weak convergence for the associated backward 
stochastic differential equation in the S-topology to derive the averaged PDE. However, 
since the averaged coefficients are discontinuous, the classical viscosity solution is not 
defined for the averaged PDE. We then use the notion of "LP— viscosity solution" in- 
troduced in [6|. We use BSDEs techniques to establish the existence of L p — viscosity 
solution for the averaged PDE. We establish weak continuity for the flow of the limit 
diffusion process and related the PDE limit to the backward stochastic differential 
equation via the representation of L p -viscosity solution. 

Keys words: Backward stochastic differential equations (BSDEs), L p -viscosity solution for 
PDEs, homogenization, S-topology, limit in Cesaro sense. 
MSC 2000 subject classifications, 60H20, 60H30, 35K60. 



1 Introduction 

Homogenization of a partial differential equation (PDE) is the process of replacing rapidly 
varying coefficients by new ones such that the solutions are close. Example: Let a be a 
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one dimensional periodic function which is uniformly elliptic. For e > 0, we consider the 
operator 

L e = div(o(j)V) 

For small e, L e can be replaced by 

L = div(aS7) 

where a is the averaged (or limit or effective) coefficient associated to a. As e — > 0, solutions 
of parabolic equations 

d t u = L £ u, u(0,x) — f(x) 

are close to the corresponding solutions with L e replaced by L. 

The probabilistic approach to homogenization gives a good description of this topic in the 
periodic or ergodic case. It is based on the asymptotic analysis of the diffusion process asso- 
ciated to the operator L e . The averaged coefficient a is then determined as a certain "mean" 
of a with respect to the invariant probability measure of the diffusion process associated to 
L. 

There is a vast literature on the homogenization of PDEs with periodic coefficients, see for 
example monographs [21 CEH EES] and the references therein. There also exists a considerable 
literature on the study of asymptotic analysis of stochastic differential equations (SDEs) with 
periodic structures and its connection with homogenization of second order partial differential 
equations (PDEs). Actually, coming from relations with semilinear and/or quasilinear PDEs 
given by a generalized Feynman-Kac formula, forward-backward SDEs (FBSDEs) have also 
been considered, see among others EH El M, M, CEH EES EDI EE] and the references therein. 

In the case where the periodicity and the ergodicity are not assumed, we don't have 
enough information about the invariant probability measure and hence, the situation is 
more delicate. 

In [13], Khasminskii & Krylov have considered the averaging of the following family of 
diffusions process 

f xY = x 1 + \Jly{x^,x^)dW s 

{ _ ( L1 ) 

{ a$> e = x 2 + j^M^M xl e )d8 + ^a^(xl' e } x 2 s ' £ )dW s 

where x\' £ is a null-recurrent fast component and xf ,e is a slow component. The function 
tp (resp. a) is IR-valued (resp. IR dx ( fc_1 ) -valued ). (W, W) is a IR fc -dimensional standard 
Brownian motion which components W is one dimensional while W is IR fc_1 -dimensional. 
They then studied the averaging of system (11.11) . They defined the averaged coefficients as 
a limit in Cesaro sense. With the additional assumption that the presumed SDE limit is 
weakly unique, they proved that the process (ex\' £ , x 2 t ,e ) converges in distribution towards 
a Markov diffusion (A^ 1 , Xj). As a byproduct, they derived the limit behavior for the linear 
PDE associated to {exl ,e , xf' £ ), in the case where the weak uniqueness holds in the Sobolev 
space W^ l loc . 

In the present note, exploiting the idea of [32], we study the homogenization of a parabolic 
semilinear PDE in the case where both the periodicity and the ergodiciy are not be assumed. 
We define the averaged coefficients as a limits in Cesaro sense. In such a way, the limit 



Homogenization of semi-linear PDEs with discontinuous coefficients 



3 



coefficients could be discontinuous. More precisely, we consider the following sequence of 
semi-linear PDEs, indexed by e > 0, 

%(s, xi, x 2 ) = £ £ (xi, x 2 )v £ (s,x 1 , x 2 ) + /(^S x 2 , v e (s, Xx, x 2 )), s G (0, t) 

(1.2) 

v £ (0, xx, x 2 ) = H(xx, x 2 ) 
where 

x d 2 ^ x d 2 x d 

C e (x l7 x 2 ) = a 00 ( — , x 2 ) — h V] Oij(— , £2)0 — 5 ^y^ & i 1} (— > ^2)^ — , 

£ <9 2 Xi z^— ' e ox 2i ox 2j £ <9x 2i 

and the real valued measurable functions / and H are defined on IR d+1 x IR and IR d+1 
respectively. 
We put, 

-if 2 := a 00 , ay := |(<r (1 VW*)y, i, j = 1, d, and a = 
One has a G ]R {d+1)xk with 

o"oo = V, 

<?oj = 0, j = 1, k - 1 
cr i0 = 0, i = 1, d 



( <p \ 
I J' 



(1.3) 



We denote, X s := (X 1 ' 6 , X 2 ' e ), b = (0, 6 (1) )*, and B = (W, W) 
The PDE (O) is then connected to the Markovian FBSDEs, 

y; = H{xt) + £ f(^f, x^, yj)d« - £ zi dM*% Vs g [o, *] 

where x = (x^x 2 ) and M xe is a martingale part of the process X s . 

It is well known that (under some conditions) the representation v e (t, x) = Y" £ holds. 

The aim of the present paper is: 

1) to show that the sequence of process (Xf,Yf,f Z^dM* e ) < s < t converges in law to 
the process (X t , Y t , f* Z u dM*) < s < t which is the unique solution to the FBSDE, 

X s = x + f°b(X u )du + f°a(X u )dB u , 0<s<t. 

(1.4) 

Y s = H(X t ) + Si f(X u , Y u )du - Jl Z u dM*, 0<s<t 

where a, b and / are respectively the average of a, b and /. 

2) As a consequence, we establish that v e tends towards v, which solves the following 
averaged equation in the //-viscosity sense. 

§j|(s, xi, x 2 ) = L(x lt x 2 )f(s,Xa, x 2 ) + J(x u x 2 , v(s, x x , x 2 )) < s < t 

(1.5) 

v(0,x x , x 2 ) = H(x l , x 2 ) 
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— d 2 - d 

where L(xu x 2 ) = / an(%i, x 2 ) — — k \~ s bAxx, x 2 ) — — , is the averaged operator. 

z — ' OXiOXj OXi 

i,j J i 

The method used to derive the averaged BSDE is based on weak convergence in the S- 
topology and is close to that used in [20]. In our framework, we show that the limit FBSDE 
( 11.41) has a unique solution. However, due to the discontinuity of the coefficients, the classical 
viscosity solution is not defined for the averaged PDE (11.51 ). We then use the notion of 
"L p — viscosity solution". We use BSDEs techniques to establish the existence of LP— viscosity 
solution for the averaged PDE. The notion of L p -viscosity solution has been introduced by 
Cafarelli et al. in [6] to study fully nonlinear PDEs with measurable coefficients. However, 
even if the notion of L p -viscosity solution is available for PDEs with merely measurable 
coefficients, one require continuity property for such solutions. In our case, the lack of L 2 - 
continuity property for the flow X x := (X 1,x , X 2,x ) transfer the difficulty to the backward 
one and hence we cannot prove the L 2 -continuity of the process Y . To overcome this difficulty, 
we establish weak continuity for the flow x 1— > {X x > x , X 2)X ) and use the fact that Y X is 
deterministic, to derive the continuity property for Yq. 

The paper is organized as follows: In section 2, we make some notations and assumptions. 
Our main results are stated in section 3. Section 4 and 5 are devoted to the proofs. 

2 Notations and assumptions 
2.1 Notations 

For a given function g(xi, x 2 ), we define 

1 f 

g + (x 2 ) := lim — / 

1 f xi 

g (x 2 ) := lim — / g(t, x 2 )dt 

Xx Jo 



g(t, x 2 )dt 



The limit in Cesaro of g is defined by, 

^{xi>o} -1- 9 {x 2 )i-{ xl <o} 



g ± (x 1 , x 2 ) := g + (x 2 )l {xi>0} +g {x 2 )h xi < 



Let p(xx, x 2 ) := a o(^i, x 2 ) 1 (= [^(p 2 (xx, x 2 )} x ) and denote by b(xx, x 2 ), a(xx, x 2 ) and 
f(xx, x 2 , y), the averaged coefficients defined as follows, 

t ( v {pb i ) ± (xx, x 2 ) . 

bi(xx, x 2 ) = — — ; — , z = l, d 

p ± {x 1 , x 2 ) 



aij{xx, x 2 ) = — I ! f 7 3 — — , i, j = 0, 1, d 



{p a ij ) ± {xx 1 x 2 ) 
p ± (x 1 , x 2 ) 

(p/)^!, x 2 , y) 



p ± {x 1 , x 2 ) 

It's worth noting that b, a and / may be discontinuous at X\ = 0. 
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2.2 Assumptions. 

We consider the following conditions, 

(Al) The function b {l \ ip are uniformly Lipschitz in the variables (x±, x 2 ), 

(A2) for each x±, their derivative in X2 up to and including second order derivatives are 
bounded continuous functions of a; 2 . 

(A3) a := (a^a^*) is uniformly elliptic, i.e: 3A > 0; V:r,£ G M d , Ca(x)£ > A|£| 2 . 
Moreover, there exist positive constants C±, C2, C3 such that 

(i) Ci < a 00 (xi, x 2 ) <C 2 

(ii) \a( Xl , x 2 )\ + |6(xi, x 2 )\ 2 < C 3 (l + \x 2 \ 2 ) 

(Bl) Let D X2 u and D 2 2 u denote respectively the gradient vector and the matrix of second 
derivatives of u with respect to X2. The following limits are uniform in X2, 



1 

xi Jo 



±1 



p(t, x 2 )dt — > p (x 2 ) as xi — > ±00 

1 f Xl 

— / D X2 p(t, x 2 )dt — > D X2 p ± {x2) as x\ — > ±00 
ci Jo 

1 f Xl 

— / Dl 2 p(t, x 2 )dt — > D 2 X2 p ± (x 2 ) as x 1 — > ±00 
^1 Jo 

(B2) For every i and j, the coefficients pfej, D X2 (pbi), D 2 X2 (pbi), pa,ij, D X2 (paij), 
D 2 2 {paij) have limits in Cesaro sense. 

(B3) For every function k G {p, ph, D X2 (pbi), D 2 X2 (pbi), pa ih D X2 (pa i:j ), D 2 X2 (pa i:j )} , there 
exists a bounded function a such that 

^ / X1 x 2 )dt - ^(xi, x 2 ) = (1 + \x 2 \ 2 )a(xi, x 2 ), 

(2.1) 



lim sup |a(xi, x 2 )\ = 0. 
x 2 el 



|xi| >00 jj^d 



(CI) 

(i) For every (x±,X2 ) G H d+1 , f(x u x 2 , •) G C 2 (IR), the bounds being uniform in (xi,x 2 ). 

(ii) if is continuous and bounded. 

(C2) pf has a limit in Cesaro sense and there exists a bounded measurable function (3 such 
that 

^ Jo" 2 p(t, aj 2 )/(t, ^2, 2/)d* - (p/) ± (^i, x 2 , y) = (1 + N 2 + |y| 2 )/3(xi, x 2 , y) 

(2.2) 

lim sup x 2 , y)\ = 0, 

|a,lHoo (x a , 1/ ) e IR <, xIR 
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(C3) For each x x , pf has derivatives up to second order in (x 2 ,y) and these derivatives are 
bounded and satisfy (C2). 

Throughout the paper, (A) stands for conditions (Al), (A2), (A3); (B) for conditions (Bl), 
(B2), (B3) and (C) for (CI), (C2), (C3). 



3 The main result. 

Consider the equation 

X f s > x = x+ [ b(X t u x )du+ I W(X% x )dB u , t<s<T (3.1) 



Assume that (A), (B) hold. Then, 

• From Khasminskii and Krylov [13] we have: the process X s := (X ,e , X 2 ' £ ) converges 
in distribution to the process X := (X 1 , X 2 ). 

and 

• From Krylov [H] we have: The limit X = (X 1 , X 2 ) is a unique weak solution to SDE 
O. 

We now define the L p -viscosity solution. This notion has been introduced by Caffarelli 
et al. in [6] to study PDEs with measurable coefficients. Wide presentation of this topic can 
be found in [12, [7]. 

Let g : [0, T] x IR d+1 x IR — > IR be a measurable function and 

Q2 _ Q 

L(xi, x 2 ) '■= > a,ij(xu x 2 ) — — h } bAxi, x 2 ) — — be the operator associated to the 

^— ' OXiOXj OXi 

i,j J i 

solution to SDE (jOl . 

We consider the parabolic equation, 

|f (s, xi, x 2 ) + L(xi, x 2 )v(s, xi, x 2 ) + g(s, x x , x 2 , v(s, x x , x 2 )) = 0, t < s <T 

(3.2) 

v(T,xi, x 2 ) = H(xt, x 2 ). 

Definition 3.1. Let p be an integer such that p > d + 2. 

(a) - A function v eC ([0, T] x IR d+1 , IR) is a L^-viscosity sub-solution of the PDE ( jOl . 
if for every x G IR d+1 , v(T, x) < H(x) and for every <p G W l p f oc ([0, T] x IR d+1 , IR) and 
(t, x) G [0, T] x IR d+1 at which v — <p has a local maximum, one has 

ess lim inf \—^-(s, x\, x 2 ) — G(s, x, (p(s, x))\ < 0. 

(s, x) { OS ) 

(b) - A function v G C ([0, T] x IR d+1 , IR) is a L p -viscosity super-solution of the PDE 
(pOl) . if for every x G IR d+1 , v(T, x) > H(x) and for every if G W^J oc ([0, T] x IR d+1 , IR) 
and (t, x) G [0, T] x IR d+1 at which v — tp has a local minimum, one has 

ess lim sup i — -?f-(s, x x , x 2 ) — G(s, x, tp(s, x)) \ > 0. 

(s,x)^(t,x) { OS ) 



Homogenization of semi-linear PDEs with discontinuous coefficients 



7 



Here, G(s, x, (p(s, x)) = L(x%, x 2 )ip(s,xi, x 2 ) + g(s, x%, x 2 , v(s, X\, X2)) is assumed to be 
merely measurable on the variable x — (x\, x 2 ). 

(c)- A function v G C ([0, T] x M d+1 , IR) is a L p -viscosity solution if it is both a W- 
viscosity sub-solution and super-solution. 

Remark 3.2. The definition (a) means that for every e > 0, r > 0, there exists a set A C 
B r (t, x) of positive measure such that, V(s, x) G A, 

dip 

— — (s, xi, x 2 ) - L(xi, x 2 )(p{s, xi, x 2 ) - g(s, x 1: x 2 , v(s, x u x 2 )) < e. 

OS 

The main result is, 

Theorem 3.3. Assume (A), (B), (C). Then, the sequence of process (X £ ,Y t £ , J Z^dM**) 
converges in law to the process (X t , Y t , f* Z u dM*) which is the unique solution to FBSDE 

O- 

Theorem 3.4. Assume (A), (B), (C). Fore > 0, letv e be the unique solution to the problem 

(CQjj. Let (Ys ) s be the unique solution of BSDE (fT^Jj. 

Then, 

(i) Equation (HE) has a unique LP -viscosity solution v such that v(t,x) = Yq\ 
(ii) for every (t,x), v £ (t,x) converges to v(t,x). 

Remark 3.5. 1) The conclusion of Theorem 13.31 remains valid if we take the forward com- 
ponent of our FBSDE (11.31) as that of [13] and replace the Brownian motion B by W. The 
dimensions and the assumptions should be accordingly rearranged. For instance, in place 
equation ( 11.31) . we consider the FBSDE, 

X! = x + j; 6(41, X^)du + ft a(*f, X^)dW u , 

(3.3) 

y; = H{xf) + Si f(^f, x* u £ , Y £ )du - Si zi dMf, Wse[o, t) 

where x = (xi,x 2 ) and 

Xl' £ = x 1 + J^(^,X^)dW s 

(3.4) 

Xl- = x 2 + J2&«(41, X^)ds + J V«(41, X^)dW s 

In this case: <p is a IR fc — valued function, W is an IR fc — Brownian motion and = (crff) 
is a d x k matrix. The nondegeneracy condition A3)-i) imposed on if should be replaced 
by, C x < £ = iV 9 «( a; i5 x 2) < C 2 . The infinitesimal generator C associated to the diffusion 
component will be more complicated since it takes account on other (mixed) second order 
derivatives. 

2) If in Theorem l3.3l we replace the initial condition of the forward component in equation 
(11.3P by (exi,x 2 ), then we obtain the same limit with initial condition (0,x 2 ) instead of 
(xi,x 2 ). 
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4.1 Tightness and convergence for the BSDE. 

Proposition 4.1. There exists a positive constant C which does not depend on e such that 

r-t 



sup IE ( sup \Y S £ \ 2 + [ \Z £ \ 2 d(M XE })\ <C. 

e \0<s<t Jo J 



Proof. Throughout this proof, K, C are positive constants which depends only on (s, t) 
and may change from line to line. It is easy to check that for all k > 1, 

sup IE ( sup [\Xl' e \ 2k + \X 2 > £ \ 2k ] j < +00. (4.1) 

£ \0<S<t J 

Using Ito's formula, we get: 

|F/| 2 + f\Zl\ 2 d{M xs ) u < \H(X*)\ 2 + K [ t \Y:\ 2 du+ f\f(xt £ ,X 2 >*,0)\ 2 du 

J s J s J s 

- 2^(Y^ZldMD- 

Passing to expectation, it is then follows by using Gronwall's lemma that, there exists a 
positive constant C which does not depend on e such that, 

ie (|r/| 2 ) < cm (\h{xd\ 2 + jf \f{x l ;\ x 2 u \ o)| 2 d^ , v* g [o, t] 

We deduce that 

IE QT \Zl\ 2 d(M x " ))j < Cm (\H{Xt)\ 2 + \f{X l ;\ X^, 0)\ 2 dn)) (4.2) 
Combining ( 14.21) and the Burkholder-Davis-Gundy inequality, we get 
Ef sup |F/| 2 + ^ [ t \Z* u \ 2 d{M xe ) u )<Cm(\H(X s t )\ 2 + f\f(gl K ,X%<, Q)\ 2 du\ 

\0<s<t * Jo J \ Jo J 

In view of condition (CI — Hi) and (14.11) . the proof is complete. ■ 

Proposition 4.2. For e > 0, let Y £ be the process defined by equation \1.3\ and M e ) be its 
martingale part. The sequence (Y £ , M 6 ) £>Q is tight on the space V ([0, t], M) x V ([0, t], M) 
endowed with the S-topology. 



Proof. Since M £ is a martingale, then by [18] or [12] . the Meyer-Zheng tightness criteria is 
fulfilled whenever 

sup ( CV(Y £ ) + m ( sup \Y £ \ + \M £ \ ) j < +00. (4.3) 

where the conditional variation CV is defined in appendix A. 
By [22], the conditional expectation CV{Y £ ) satisfies 

•t 



cv(y £ ) < Km n \f(x]> £ , x 2 ' £ , r/)|W) , 



where K is a constant which only depends on t. 

Combining condition (CI) and Proposition 14.11 we deduce (14.31) . 
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Proposition 4.3. There exists (Y, M) and a countable subset D of [0, t] such that along a 
subsequence of e, 

(i) (Y £ ,M £ ) ^ (Y, M) onV([0, t], M) xD([0, t], M) endowed with the S-topology. 

(ii) (Y e , M £ ) — > (Y, M) in finite- distribution on D c . 

(Hi) (X 1,£ ,X 2 ' £ ,Y £ ) => (X 1 , X 2 ,Y) , in the sense of weak convergence in c(m+ x m d+1 ) x 
D(M + x IR), equipped with the product of the uniform convergence on compact sets and the 
S topologies. 

Proof, (i)- From Proposition321 the family (y e , M £ ) £ is tight on V ([0, t], R)xP([0, t], M) 
endowed with the S-topology. Hence, along a subsequence (still denoted by e), (Y £ , M £ ) £ 
converges in law on V ([0, t], IR) x D ([0, t], IR) towards a cad-lag process (Y, M). 

(ii)- Thanks to Theorem 3.1 in Jakubowski [12], there exists a countable set D such 
that along a subsequence the convergence in law holds. Moreover, the convergence in finite- 
distribution holds on D c . ■ 



4.2 Identification of the limit finite variation process. 

Proposition 4.4. Let (Y,M), the limit process defined in Proposition ^. 3\ Then, 

(i) For every s G [0, t] — D, 

Y s = H(X t ) + fJ(Xl XI Y)du - (M t - M s ), 

(4.4) 

iE(sup < s < t |r s | 2 + |x|| 2 + |x s 2 | 2 ) <c. 

(ii) M is a T s -martingale, where T s '■= cr{X u , Y u , < u < s} augmented with the IP-null 
sets. 

To prove this proposition, we need the following lemmas. 

Lemma 4.5. Assume (A), (B), (C2) and (C3). For y G IR, let V £ (xi 1 x 2 ,y) denote the 
solution of the following equation: 

o o(t' x 2)D 2 Xl u(x u x 2 ,y) = f(f, x 2 , y) - f(x h x 2 , y) 

(4.5) 

u(0, x 2 ) = D Xl u(0, x 2 ) = 0. 
Then, 

(i) D Xl V £ (x 1 , x 2 ,y) = xi(l + \x 2 \ 2 + \y\ 2 )(3(— , x 2 , y) - xi(l + \x 2 \ 2 )m( — ,x 2 , y), 

(ii) foranyK £ (x u x 2 ,y) G {V £ , D X2 V £ , D 2 X2 V £ , D Xl D X2 V £ , D y V £ , D 2 V £ , D Xl D y V £ , D X2 D y V £ } 
it holds, 

X\ X\ 

K £ (x u x 2 ,y) = x\(l + \x 2 \ 2 + \y\ 2 )p( — , x 2 , y) + x\(l + |x 2 | 2 )m( — , x 2 , y) 

where m(^r, x 2 , y) := - ^ [ x l] x l) ~ a ( i~ ' X2 ) an< ^ ^2), x 2 , y) are various bounded 

functions which satisfy property \2. l\) and \2.2\) respectively. 
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Proof. We will adapt the idea of [13] to our situation. For a fixed y, we set 

F(— , x 2 ,y) ■= — I p(~, x 2 )g(-, x 2 , y)dt 
e x\ Jo e e 

where g{f, x 2 , y) := f(f, x 2 , y) -J(xi, x 2 , y). 
For X\ > 0, we have 



F(—, x 2 ,y) 



1 f Xl t t 

— / p(-, x 2 )f(-, x 2 , y)dt- (pf) + (x 2 , y) 
xi Jo £ £ 

(pf)+(x 2 , y) 1 



+ (pf) + (x 2 ,y) 



p + (x 2 ) Xi Jo £ 



p(-, x 2 )dt 



,X\ 



;i + |x 2 | 2 + br)A(— , x 2 , y ) 



+ (Pf) + (x 2 ,y) 



1 - 



P + (x 2 )xi Jo e 



p(-, x 2 )dt 



,Xi 



[1 + \x 2 \ + \y\ )Pi(—, x 2 , y)-(l+ \x 2 \ 



(pf) + (x 2 , y) x 1 

— - — a x {— , x 2 ) 

P + {x 2 ) £ 



Since, D Xl V e (xi, x 2 ,y) = x\F(2^, x 2 ,y), we derive the result for D Xl V e (xi, x 2 ,y). Further, 
by integrating, we get 



V £ {xi, x 2 ,y) 



x\(l + \x 2 \ 2 + \y\ 



Xl Jo 



'1 + \x 2 \ 



(pf) + (x 2 , y) e 



p + (x 2 ) 



X\ Jo 



t/9i(t, x 2 , y)dt 
x 2 )dt 



Clearly, (3(f, x 2 , y) := (^) 2 / £ W*, x 2 , y)dt, a(f, x 2 ) = f~< ta(t, x 2 )dt 
satisfy (12.10 and (12.21) respectively. The same argument can be used in the case X\ < 0. The 
result for D X2 V e (x 1 , x 2 ,y), D 2 2 V £ (xi,x 2 ,y) and D Xl D X2 V e (xi, x 2 , y) can obtained by using 
similar arguments. ■ 



Lemma 4.6. sup 

0<s<t 

in probability as e — > 0. 



/(^!, X 2 U > £ , Y £ ) - fW, X 2 U > £ , Y £ ) ) du 



tends to zero 



Proof. Let V s denote the solution of equation (14.51) . Note that V s has first and second 
derivatives in (xi, x 2 ,y) which are possibly discontinuous only at X\ = 0. Then, as in [13], 
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since if is non degenerate, we can use Ito-Krylov formula to get 



X 



he 



V{X?, X 2 ' £ , Y/) = V% Xl , x 2 ,Y £ ) + / [/(-*-, X 2 ' £ , FJ) - /(X^, X 2 > £ , 



<9 2 F £ 



dx 2i dx 2j 



■Mfvh 



-2,e 



dV e 



dx 



2.1 



/ s p>ye aye , 

i^(xl' £ , x 2 u £ , y^mx^, x^)dw u + ^(xl' £ , x 2 u > £ , n> (1) PC' E , X^)dW u \ 



1 

+ - 



dV e X 1,£ 

^pe, x 2 /, van-*-, x 2 /, r:)du + 

o °y 6 

s pflye Y 1 ' 6 



dV £ X^-' £ 
-^(X^,X^,Y^)Z £ a(^,X^)dB u 

o °y 6 



2 Jo dydy 



Xu' £ )(Z £ u Tdu 



(4.6) 



+ 



2 Jo dydx ' ~ : 

In view of Lemma [4751 it is obvious to see that V e (xi, x 2 ,Y Q £ ) tends to zero. Once again, 
from Lemma [4.51 we have 



\V%Xl>*,X^,Y*)\ 



< 



1,6 



(1 + |XM 2 +|F/| 2 )|/5(^,X 2 ' £ ,F/ 



+ e 



(1 + |X 2 ' £ | 2 

{|Xi' E |>v^} 



( P /) ± (x 2 . £ ,y/) l _ / x s M ^ 



a 



(l + |X 2 ' e | 2 + |F/| 2 )|/?(^ ^ 2 - £ ^- 



+ 1 {\xl' s \>VE}\ X s 



l,e|2 



(1 + IX 



2,e|2 



(p/m s 2 ' £ ,n £ ),„^. 



l,e 



2,£\ 



2,e\ 



From the estimates of the processes X]' £ , X 2 ' £ , F/ and the fact that (p/) 1 * 1 satisfies conditions 
(C), we deduce that 



0<s<t 



X 



,x 



IE sup \V £ ' n {Xl' £ , X 2 > £ ,y/)| < K \e+ sup sup \/3(—, x\ y)\ + sup sup|a( — , x 



|xi|>Ve (x 2 ,y) 



Then, since a and (3 satisfy (12. ip and (12. 2p respectively, the right hand side of the previous 
inequality tends to zero as e — > 0. Similarly, one can show that 



d 2 V £ 



+ / aij(X u ' £ , X 2 '^ 



(X^ £ ,Xr,Y:)du+ I bf\x l u > £ ,Xl 



-2,E 



dV £ 



dx 



2.7 



{X^,X^,Y £ )du 



dx 2i dx 2j 

/ s 8V £ dV £ ~ 

[^-(X^' £ , X 2 < £ , Y £ MX^, X^)dW u + -^(X 1 ^, X 2 ' £ , Y^\X^, X 2 u £ )dW u ] 



dV £ X 1,£ 

-dy-( X u £ i X u £ i Y:)f(^ X 2 * £ , Y £ )du + 



dV £ X 1 ' 6 
■ -i-(Xl' £ ,X^,Y:)Z:a(^,Xl' £ )dB u 

o °y e 



i 

+ ■x 



s pfiye Yhe 
u v '"l.e v 2,e ve\r 7 e^*( yv u 



2 J dydy 



(X u ' £ ,X u ' e } Y £ )Z^aa* 



x 2 /){z £ u ydu 



1 

+ - 



^ \x^,X 2 ' £ ,Y:)aa*(^,X 2 n(Z £ u ydu 



2 Jo dydx 
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converge to zero in probability. Let us give an explanation concerning the one but last term, 
which is the most delicate one. 



o dydy 



<C sup (1 + \X 2 U £ \ 2 ) 



0<u<s 



d 2 V £ 



dydy 



( yl,e y2,e v e\ 

l A « 5 A M 1 Y U, 



\Zi\ 2 du 



bince {/„ \Z £ u \ 2 du, < s < t} is the increasing process associated to a bounded martingale, 
so the L P (IP) norm of J * \Z^\ 2 du is bounded, for all p > 1. Moreover sup 0<u<t (l + |X 2,£ | 2 ) is 
also bounded in all L P (JP) spaces. Finally the same argument as above shows that 

d 2 V £ 



sup 

0<u<s 



dydy 



as e 



in probability, as e — > 0. 

Lemma 4.7. / ' J(X^ , X 2 U > £ , Y £ )du ^ [ 'j(Xl X 2 , Y u )du on C([0, t], M) 
Jo Jo 

For the proof of this Lemma, we need the following two results. 

Lemma 4.8. Let X] := x x + \ ^>{X\, X 2 )dW u , < s < t, and, assume (A2-i), (Bl). 

Jo 

Fore> 0, let D £ n : = js : s G [0, t] / X 1 / G 5(0, ^)|. 

Define also D n :=|s:se [0, t] / X] G B(0, ^)|. 

Then, there exists a constant c > such that for each n > 1, e > 0, 



0. 



IE\D £ n \ < 



and 



IE\DJ < 



n 



n 



where |. | denotes the Lebesgue measure. 

Proof. Consider the sequence ($? n ) of functions defined as follows, 

—x/n—l/2n 2 if x<—l/n 

V n (x) = I x 2 /2 if -l/n<x <l/n 
x/n—\/2n 2 if x>l/n 

We put, Tp := a 00 := p(x u x 2 y 1 . 
Using Ito's formula, we get 

9 n (Xl) = 9 n (X^) + J* ^' n (Xl)Tp(Xl, X 2 )dW s + ± jf *" n (X\)7?(Xl X 2 )ds, s G [0, t] 
Since Tp is lower bounded by Ci, taking the expectation, we get 

dIE f \ { _,_ M {X\)ds <W> f\l{Xl)T5 2 (XlX 2 s )d S 
Jo n n Jo 

= 21E [* n (XD - * n ( Xl )] 

It follows that E(|D n |) < 2C 1 ~ 1 IE [^(X*) - * n (xi)] < c/n. The same argument, applies to 
-D^, allows us to show the first estimate. ■ 
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Lemma 4.9. Consider a collection {Z £ , e > 0} of real valued random variables, and a real 
valued random variable Z. Assume that for each n> 1, we have the decompositions 

z £ = z^ + zft 

Z — Z n + Zt.i 



such that for each fixed n > 1, 



zl 



m\z 2 f\< ' 



m\zi\ < 



n 
c 



n 



Then Z £ Z , as e -> 0. 

Proof. The above assumptions imply that the collection of random variables {Z e , e > 0} 
is tight. Hence the result will follow from the fact that 

E$(Z £ ) -> IE$(Z), as e -> 

for all $ G Ct,(IR) which is uniformly Lipschitz. Let $ be such a function, and denote by K 
its Lipschitz constant. Then 

|E$(Z £ ) - m$(Z)\ < E|$(Z £ ) - $(^- £ )| + +|E<I>(^' £ ) - + - $(Z)| 

< |E$(Z^ £ ) - + 2AT- " 



jn 
Hence 

limsup |IE$(Z £ ) - < 2^^, 

for all n > 1. The result follows. ■ 

Proof of Lemma 14. 7L For each n > 1, define a function 6 n G C(H, [0, 1]) such that 
9 n {x) = for \x\ < l/(2n), and # n (a;) = 1 for |x| > 1/n. We have 

7(Xl' £ , X s 2 ' £ , Y!)ds = f f{X?, X 2 /, Y!)9 n {X x /)ds + f J(X^, X s 2 < £ , Y!) [1 - 9 n (X^)]ds 

Jo Jo 

— + 

J(Xl,X 2 ,Y s )ds = fj{XlXlY s )6 n {Xl)ds + [ J(Xl,X',Y s )[l - e n {X])]ds 
Jo Jo 

= zl + zl 

Note that the mapping 

(x\x 2 ,y)^ / J(xl,x 2 s ,y s )9 n (xl)ds 



is continuous from C([0, t]) x D([0, t]) equipped with the product of the sup-norm and the 
S topologies into IR. Hence, from Proposition 14.31 Z^ 6 =>- Z\ as e — > 0, for each fixed n > 1. 
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Moreover, from Lemma l4~8l the linear growth property of /, Proposition 14.11 and (14.11) . we 
deduce that 



E\Z 2 f\ < -±=, E\Z 2 n \ < 



Lemma [4.71 now follows from Lemma \4 .91 ■ 

Proof of Proposition 14.41 Passing to the limit in the backward component of the equation 
A1.3B and using Lemmas 14.61 and 14.71 we derive assertion (z) . 

Assertion (ii) can be proved by using the same arguments that developed in [2T] section 6. 



4.3 Identification of the limit martingale. 

Proposition 4.10. Let (Y S ,Z S , < s < t) be the unique solution to BSDE (H(X t ), /). 
Then, for every s G [0, t], 

IE\Y S -Y s \ 2 + lE(lM- jf Z u dM*] t - [M - jf Z u dM*] s j = 0. 

Proof. For every s G [0, t] — D, we have 

Y s = H(X t ) + fj(X u , Y u )du - (M t - M s ) 

Y s = H(X t ) + fJ(X u , Y u )du - /* Z u dM* 

Arguing as in [2T], we show that M :— f s Z u dM^ is a jF s -martingale. 
Since / satisfies condition (CI), we get by Ito's formula, that 

m\Y a - Y s \ 2 + IE {[M - jf Z u dM* ] t - [M - jf Z tt dM* < C1E £ \Y U - Y u \ 2 du. 

Therefore, Gronwall's lemma yields that IE|F S — Y s \ 2 = 0, Vs G [0, t] — D. But, F is 
continuous, F is cad-lag, and, D is countable. Hence, Y s = Y s , W-a.s, Vs G [0, t]. 

Moreover, we deduce that, IE \ [M - Z u dM*\ t - [M - J Z u dM*]^J =0. ■ 
As a consequence of Proposition 14.101 we have 

Corollary 4.11. (y £ , jf Z*dM**J ^ (y, Z u dMfj. 



5 Proof of Theorem I3.4L 



Since, under assumptions (Al) and (A2), the SDE (|3.1I) is weakly unique, the martingale 
problem associated to X = (X 1 , X 2 ) is well posed. We then have the following: 

Proposition 5.1. Assume that g satisfies (CI). Then, 
(i) For a fixed positive number T, the BSDE 

rt, x i^irX 1 ** 



Y*>* = H(X%*)+ f g(u,X^,Y^)du- f ft 

J s J s 



dM* ' ,t<s<T. 
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admits a unique solution (Y* ,x , Zl> x ) t < s <T such that the component (Yg ,x )t< s <T is bounded. 

{ii) If moreover, the function (t, x) G [0, T] x M d+1 i — ► v(t, x) := Y^' x is continuous, then 
it is a LP -viscosity solution of the PDE jl3.S\) . 

Proof, (i) Thanks to Remark 3.5 of [20], it is enough to prove existence and uniqueness for 
the BSDE, 

Y?* = H{X%*) + J T g(u, Xi', Y^)du - £ Z^dB u ,t< s < T. 

But, this can be proved by usual arguments of BSDEs. For instance, it's obvious that 
uniqueness holds under (CI), and, we can prove the existence of the solution by using a 
Picard type approximation. 

(ii) We only prove that v is L p - viscosity sub-solution. The proof for super-solution can be 
performed similarly. For <p G W^ oc ([0, T] x IR d+1 , Hi), let (t, x) G [0, T] x M d+1 be a point 
which is a local maximum of v — ip. Since p > d + 2, then ip has a continuous version which 
we consider bellow. We assume without loss of generality that v(t, x) = ip(t, x). Following 
Cafarelli et al., assume that there exists e , r > such that 

§£(s, xi, x 2 ) + L(xi, x 2 )ip(s,x 1 , x 2 ) + g(s, xi, x 2 , v(s, x 1: x 2 )) < -e , X-a.s in B ro (t, x) 
v(s, x) < (p(s, x) + e (s — t), X-a.s in B ro (t, x) 
Let Aq G B ro (t, x) be a set of positive measure such that (t, x) G A . Define 

r = inf {s > t; (s, Xl' x ) £A }A (t + r ) 

The process (F s , Z s ) = ((Y^%), l [0 ,r](s)(Z^')) se [F,f+n,] solves then the BSDE 

Y s = Vl {r, X* T >*) + / l [0 , T ](u)g(u, X^, v{u, X$*))du - / Z u dM?", s G [t, t + r }. 

J s J s 



On other hand, setting ip(s, x) = cp(s, x) + e (s — t), we have by Ito-Krylov's formula that 

^ sefct+ro] 

»t+r a nt+r 



the process (Y s , Z s ) = (^(s, X^), l [0jT] (s)V^(s, X^)) t ~ solves the BSDE 



Y s = ^(r, Xl>*) - J l [0 ,r](«)[eo + (^ + L<p)(u, X^))du - J Z u dM u 



>From the choice of r, (r, X*' x ) G A . Therefore v(r, X^< x ) < X^) and thanks to the 
comparison theorem [20], we deduce that Y^ < Yj?, i.e vi(t, x) < <p(t, x), which contradicts 
our hypothesis. ■ 

Remark 5.2. (i) Whenever g does not depends on t; v(t, x) = Y(f is a //-viscosity solution 
of the PDE 

§jf(s, xi, x 2 ) = L(x u x 2 )v{s, xx, x 2 ) +g(x lt x 2 , v(s, x u x 2 )) 



v(0, xi, x 2 ) = H(x±, x 2 ), s > 0, x = (x±, x 2 ) G IR' 



d+l 
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where (X x , Yf,Z x ; < s < t), solves the following decoupled FBSDE 
X* = x + ft b(X*)du + j; a(X x u )dB u , 0<s<t. 

Y* = H(Xf) + f s g{XZ, Y u x )du - Z x dM^ , 0<s<t 

(ii) Since / satisfies (C) and p is bounded, one can easily verify that / satisfies (CI). 
Therefore, for a fixed positive t, the BSDE with data (H(Xf), f) admit a unique solution 
(Yf, Z x ) < s <t- Moreover, if the function (i, x) G IR+ x IR d+1 i— > v(t, x) = Yq is continuous 
and hence, it is a L p - viscosity solution of the PDE (11.51) . 

Under assumptions (A), (B), the SDE ( 13. ip has a unique weak solution [IB]. We then 
have, 

Proposition 5.3. (Continuity in law of the flow x \— > X x ) 

Assume (A), (B). Let X* be the unique weak solution of the SDE (\3.1\) . and 

X:-.= x n + [ b(X™)du+ [ a(X%)dB u , 0<s<t 
Jo Jo 

Assume that x n converges towards x = (x 1 , x 2 ) G IR 1+d . Then, X n X x . 

Proof. Since b and a satisfy (A), (B), one can easily check that the sequence X n is tight 
in C([0, t] x IR d+1 ). By Prokhorov's theorem, there exists a subsequence (denoted also by 
X n ) which converges weakly to a process X. We shall show that X is a weak solution of 

sde dam . 

• Step 1: For every <p G C™(SR 1+d ), 

Vw G [0, t], <f(X u ) — / Lip(X v )dv is a jF x -martingale. 

Jo 

We have need to show that for every <p G C£° (IR 1+a! ) , every < s < u and every function § s 
of (X Xn ) < r<s which is bounded and continuous in the topology of the uniform convergence, 

= IE jfoW) - ¥>W") - £l<p(X?)dv]$.(X*"] 

-±>m^[<p(x u )-<p(x a )- Lip(x v )dv]$ a (x. 

Indeed, since <p, $ are continuous functions and L is continuous out of the set {x\ = 0}, 
similar argument as that developed in the proof of Lemma T4.7I gives 

[<p(X?) - <p{X?) - / T^{X?)dv]*.{X*~) ^ [p{X u ) - <p(X s ) - / Lcp(X v )dv}<S> s (X.) 



Since cp, $ are bounded functions and sup n IE(sup sg [ t i \X Xn \ ) < oo, the result follows by the 
uniform integrability criterium. Hence, Ie|[<^(X u ) — <p(X s ) — J L<p(X v )dv)<& s (X .) | = 

and therefore <p(X u ) — p(X s ) — f u Lp(X v )dv is a jF x -martingale. 

• Step 2: From step 1, there exists a jF x -Brownian motion B such that, 

X s = x + [ b(X u )du + / a(X u )dB u , < s < t. 
Jo Jo 



Weak uniqueness for SDE (OLE]) allows us to deduce that X x " ^ X x . 
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Proposition 5.4. Assume (A), (B), (C). Then, 
(<) limF* =F ^. 

(m) T/je map (t,x) i — > F t,:r is continuous. 



(Hi) For p > d + 2, the function v(t,x) := F ,:c is a L p -viscosity solution to the PDE ( ti.51) . 

Proof, (i) Let F be the limit process defined in Proposition (14.31) . Since F e and F are 
deterministic, it is enough to prove that lim £ ^ E(F e ) = IE(F ). We have, 

F E = H{Xf) + Ji f(T u , X^, Yt)du - Ml 

F = H(X t ) + f*J(X U7 Y u )du - M t 

>From Jakubowski [12], the projection: y i— > y t is continuous in the S-topology. We then 
deduce that F e converges towards F in distribution. Since F £ and F are bounded, then 
hm £ ^ IE(F £ ) =E(Fo). 

(ii) Let (t n , x„) — > (£, ar). We assume that t > t n > 0. We have, 

Y t n ,* n = H (X£) + f(X^,Y^)du- Zi^dMf\ 0< s< tn, (5.1) 

J s J s 

where X Xn ^ X*. 

Since if is a bounded continuous function and / satisfies (CI), one can easily show that the 
sequence {(Y*-*», / Q l [Sjtn] («)^dM^ n )} neIN . is tight in V([0, tjxMx JR). 
Let us rewrite the equation (15.11) as follows, 

yt n , Xn = H (X£) + J f(X^,Y^)du- J* l [sM {u)Zi^dMr n (5.2) 

J(X^,Y^)du,0< s<t. 
At + A 2 n 



Convergence of A 



2 



One has IE 



f{X?, Y^ Xn )du 



< K{\x\)\t — t n \. Hence A 2 n tends to zero in probability. 



• Convergence of A n 

Denote by (F', M') the weak limit of {(Y tn > Xn , J Q \ [sM (u)Z^dM^ Xn )} n&w . In view of 



Lemma[47l one has I f(X* n , Y*"' Xn )du ^ I f(X*, Y^)du. 

J S J S 

Passing to the limit in (15.21) . we obtain that 

f; = H(xf ) + J J(x*, Y^du - (m; - m' s ), 8 e [o, t] n d c . 



The uniqueness of the considered BSDE ensures that Vs G [0, t], Y' a = F s t,!E P-ps. Hence 
Y tn ' Xn ^ Y t,x . As in (i), one derive th; 

v t,x 

Assertion (Hi) follows from Remark I5.2[ 



ytn,x n yt,X_ Ag ^ deriye that yt n ,x n l™ yt, X ^ ^ continuity Q f 
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Remark As in KK, we can take W instead of W . 



A Appendix: S-topology 

The S-topology has been introduced by Jakubowski ([12], 1997) as a topology defined on 
the Skorohod space of cadlag functions: T>([0, T]; IR). This topology is weaker than the 
Skorohod topology but tightness criteria are easier to establish. These criteria are the same 
as the one used in Meyer-Zheng topology, ([18], 1984). 

Let N a,b (z) denotes the number of up-crossing of the function z G ^([0, T]; IR) in a given 
level a < b. We recall some facts about the S-topology. 

Proposition A.l. (A criteria for S-tight). A sequence (Y £ ) e>0 is S-tight if and only if it is 
relatively compact on the S-topology. 

Let (Y £ ) £>0 be a family of stochastic processes in T>([0, T}; IR). Then this family is tight for 
the S-topology if and only if (||^ £ ||oo)e>o an d {N a,b (Y £ )) e>0 are tight for each a < b. 

Let (ft, T, IP, (T t )t>o) be a stochastic basis. If (Y)o<t<T is a process in T>([0, T}; IR) such 
that Y t is integrable for any t, the conditional variation of Y is defined by 

n-1 

CV(Y) = sup J]lE[|IE[y ti+1 - Y u | F ti ]\]. 

0<ti<...<t n =T, partition of [0,T] . 1 

The process is call quasimartingale if CV(Y) < +oo. When Y is a JF r martingale, CV(Y) = 
0. A variation of Doob inequality (cf. lemma 3, p. 359 in Meyer and Zheng, 1984, where it 
is assumed that Yt = 0) implies that 



IP 


sup \Y t \ > k 


< | (cV(Y) + m 


sup \Y t \ 




*e[o,T] 




te[o,T] 



IE [N a ' b (Y)] < 



\a\ +CV(Y) + IE 



sup | F t | 
te[o,T] 



It follows that a sequence (Y £ ) £>0 is S-tight if 

CV(Y £ ) + IE 



sup 

£>0 



sup \Y t e 
te[o,T] 



< 



-oo. 



Theorem A. 2. Let (Y £ ) £>0 be a S-tight family of stochastic process in V([0, T}; IR). Then 
there exists a sequence (^fc) fcg j/v decreasing to zero, some process Y G T>([0, T\\ IR) and a 
countable subset D G [0, T] such that for any n and any (t\, t n ) G [0, T]\D, 



Y 

•••5 1 t 



) ^ (Y tl , Y tn ) 



Remark A. 3. The projection :tc t y G (T>([0, T]; IR), S) i— > y(T)is continuous (see Remark 
2.4, p. 8 in Jakubowski, 1997), but y h- > y(t) is not continuous for each < t < T. 

Lemma A. 4. Let (U £ , M £ ) be a multidimensional process in ^([0, T]; IR P ) (p G IN*) con- 



verging to (U, M) in the S-topology. Let (T\ 



)t>o 



(resp. (J-f )t>o) be the minimal complete 



admissible filtration generated by U £ (resp. U ). We assume moreover that, for every T > 0, 
sup £>0 IE [sup < t < T |Mf| 2 ] < C T . 

If M £ is a J- -martingale and M is T v -adapted, then M is a T v -martingale. 



Homogenization of semi-linear PDEs with discontinuous coefficients 



19 



Lemma A. 5. Let (Y £ ) e>0 be a sequence of process converging weakly in T>([0, T]; M p ) to Y . 
We assume that sup e>0 IE [sup 0<t<T |l^ e | 2 ] < +00. Hence, for anyt > 0, E [sup 0<i<T \Y t \ 2 ~\ < 
+00. 
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